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ABSTRACT 

We review some important topics related to the semirelativistic description of bound 
states by the spinless Salpeter equation: the special case of the Coulomb interaction, 
numerical approximation methods, and a way to avoid the problematic square- root 
operator of the relativistic kinetic energy. 



1. Introduction 



The "spinless Salpeter equation" represents a well-defined standard approximation 
to the Bethe-Salpeter formalism for the description of bound states within relativistic 
quantum field theories. It may be derived from the Bethe-Salpeter equationB 

1. by eliminating, in full accordance with the spirit of an instantaneous interaction, 
any dependence on timelike variables, which leads to the Salpeter equation,B and 

2. by neglecting any reference to the spin degrees of freedom of the involved bound- 
state constituents and restricting exclusively to positive-energy solutions. 

The Hamiltonian governing the dynamics of the quantum system under consideration 

• incorporates relativistic kinematics by involving the square-root operator of the 
relativistic kinetic energy, \/p'^ + m? for particles of mass m and momentum p, 
but 

• describes the forces acting between the bound-state constituents by an arbitrary 
coordinate-dependent static interaction potential ^(x). 

For the case of bound states consisting of two particles of equal mass m, the generic 
Hamiltonian H in the center-of-momentum frame of these constituents, expressed in 
terms of the relative momentum p and the relative coordinate x, reads 



if = 2Y'p2 + m2 + 1/(x) . (1) 
The equation of motion involving this type of Hamiltonian, 

H\ij) = E\i;) , 

with energy eigenvalues E and corresponding Hilbert-space eigenstates \ip), has been 



widely used, e. g., for the (semi-)relativistic descripj 
of quarks within the framework of potential models.^ 



of hadrons as bound states 



2. The Spinless Relativistic Coulomb Problem 

Without doubt, a central role in physics is played by the Coulomb potential Vc(r). 
This Coulomb potential is a spherically symmetric potential, i. e., one which depends 
only on the radial coordinate r = |x|; its interaction strength is parametrized by some 
coupling constant n: 

y(x) = rc(r) = -- , «:>0. (2) 
r 

The bound-state problem defined by the semirelativistic Hamiltonian of Eq. (|I|) with 
the Coulomb potential (|^) is what we call the "spinless relativistic Coulomb problem." 



2.1. (Mostly) Analytic Results 

Over the past years, the spinless relativistic Coulomb problem has been the subject 
of intense study. We summarize, in chronological order, the knowledge gained so far: 
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Herbst,Q in a rigorous mathematical discussion, developed the complete spectral 
theory of the one-particle counterpart of the operator ([1|), (Q) — from which one 
may directly deduce for the two-particle semirelativistic Coulombic Hamiltonian 
under consideration its essential self-adjointness for k < 1 and the existence of 
its Friedrichs extension up to some critical value, viz., 

_ 4 

of the coupling constant — and derived some strict lower bound on the energy 
Eq of the ground state which translates in the two-particle case to 



Durand and DurandS presented in an involved analysis of the spinless relativistic 
Coulomb problem for the particular case of vanishing orbital angular momentum 
of the bound-state constituents the explicit construction of the analytic solution 
for the corresponding wave function. 

Castorina et a/.i generalized the behaviour of this wave function near the origin, 
that is, for small relative distances of the bound-state constituents, to arbitrary 
values of the orbital angular momentum. 

Hardekopf and SucheJl carried out a comprehensive numerical investigation of 
one- and two-particle relativistic wave equations for both spin-0 as well as spin-^ 
particles. Their numerical results for the ground-state energy Eq of the spinless 
relativistic Coulomb problem are consistent with the lower bound of Herbst and 
provide no evidence that this ground-state energy indeed approaches zero as the 
coupling constant k rises to its critical value = 4/77, that is, 

Eo{k = Kcr) 7^ . 

Martin and RoyS improved the lower bound of Herbst somewhat to 




Raynal et al^ made use of a generalized version of the "local-energy" theorem:lli 
Assume 

1. that the Fourier transform V^(p) of the interaction potential V^(x) is strictly 
negative, except at infinity, as is the case for the Coulomb potential (|), 

2. that the spectrum of the Hamiltonian under consideration, H, is discrete, 
and 
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3. that the ground state of the Hamiltonian exists. 



Define the local energy 

/c/3g^(p_q)0(q) 



S{p) ^ T(p) + 



IP) 



here T(p) represents the kinetic energy, which in our case is, of course, given 
by 



T(p) = 2^p2 + m2 , 
and 0(p) denotes some suitably chosen, positive trial function, 

0(P) > . 

Then the lowest-lying eigenvalue of the Hamiltonian H, Eo, is bounded by 

inf£^(p) < Eq < sup£^(p) . 
P p 

With the help of this theorem, Raynal et al. succeeded in restricting numerically 
the ground-state energy eigenvalue of the semirelativistic Hamiltonian (|l]), (H), 
considered as a function of the coupling strength k, to some remarkably narrow 
band. In particular, at the critical value = 4/7r of this coupling constant n 
they found 

0.4825 < ^"^"^ = "^"^ < 0.4843 , 
2 m 

which proved the nonvanishing of the ground-state energy eigenvalue Eq at the 
critical coupling constant. 

Le Yaouanc et a/.0 gave the systematic series expansion of the eigenvalues of the 
Hamiltonian (|l]), (j^) in powers of the coupling constant k up to and including 
the order 0{k^) for arbitrary states with vanishing orbital angular momentum. 
For the ground-state energy Eq this expansion reads 

Eo ^ tc' k\ ^. g. 



2m 8 128 12 vr 64 

Higher orders of this expansion in the form given by Le Yaouanc et al., however, 
involve explicitly the corresponding wave functions and their derivatives, which 
renders difficult to derive any profit from it. Up to order 0{k,^) the above result 
coincides with the improved lower bound of Martin and Roy. 

The presence of non-analytic In k terms in the expansion quoted by Le Yaouanc et al. 
gives a clear hint at the non-analytic nature to be expected for the energy eigenvalues 
of the spinless relativistic Coulomb problem as functions of the coupling strength k. 
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2.2. Analytic Upper Bounds 

In view of the fact that exact solutions to the spinless relativistic Coulomb problem 
are still lacking, we seek for an analytic upper bound on the ground-state energy level 
of the semirelativistic Coulombic Hamiltonian (|l|), (0). To this end, we make use of a 
rather standard variational technique, which proceeds along the steps of the following, 
extremely simple recipe:@ 

1. Choose a suitable set of trial states {|A)}. The different members of this set are 
distinguished from each other by some sort of variational parameter A. 

2. Compute the set of expectation values of the Hamiltonian under consideration, 
H, with respect to these trial states |A) in order to obtain 

E{X) ^ {X\H\X) . 



3. Determine, from the first derivative, that value Amin of the variational parameter 
A which minimizes the resulting, A-dependent expression E{X). 

4. Compute E{X) at the point of the minimum Amin to find in this way the minimal 
expectation value E{X^i^) of the Hamiltonian H in the Hilbert-space subsector 
of the chosen trial states |A). 

This minimum ii^(Amin) provides, of course, only an upper bound to the proper energy 
eigenvalue E of the Hamiltonian H: 

E < E{Xinm) ■ 

Let us apply the above simple recipe to the case of the Coulomb potential, Eq. (H). 
For the Coulomb potential, the most reasonable choice of trial states is obviously the 
one for which the coordinate-space representation tp{x) of the states |A) for vanishing 
radial and orbital angular momentum quantum numbers is given by the hydrogen-like 
trial functions 

V'(x) = W — exp(-Ar) , A > . 

For this particular set of trial functions we obtain for the expectation values we shall 
be interested in, namely, the ones of the square of the momentum p and of the inverse 
of the radial coordinate r, respectively, evaluated with respect to the trial states |A): 

(A|p2|A) = A^ 

and 

' a\ = A 



A 

Now, as an immediate consequence of the fundamental postulates of any quantum 
theory, the expectation value of some given Hamiltonian H taken with respect to any 
normalized Hilbert-space state and, therefore, in particular, taken with respect to any 
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of the above trial states must necessarily be larger than or equal to that eigenvalue 
Eq of the Hamiltonian H which corresponds to its ground state: 



Eo < E{\) = {\\H\\) . 

The application to the semirelativistic Hamiltonian of Eq. (|l]) yields for the right-hand 
side of this inequality 

E{X) = 2 ^p^ + m2 + (A|V^(x)|A) . (3) 

In order to obtain a first crude estimate, we may take advantage of some trivial but 
nevertheless fundamental inequality. This inequality relates the expectation values of 
both the first and second powers of any self-adjoint but otherwise arbitrary operator 
= 0^ taken with respect to arbitrary Hilbert-space vectors |) (in the domain of O) 
normalized to unity; it reads 

\m<^). 

Applying this inequality to the kinetic-energy part of the above expression for E{X), 
we may replace, in turn, E{X) by an upper bound which can be evaluated much easier 
than E{X) itself: 



^(A) < 2^(A|p2|A) + m2 + (A|V(x)|A) . 

Identifying in this — as far as its evaluation is concerned, simplified — upper bound the 
until-now general potential V{x.) with the Coulomb potential (0) and inserting both 
of the A-dependent expectation values given above implies 



E{X) < 2VX^ + - kX . 

From this intermediate result, upon inspection of the limit A — 00, we may state 
already at this rather early stage that, for the semirelativistic Hamiltonian (|^) to 
be bounded from below at all, the Coulombic coupling strength k has to stay below 
a certain critical value: 

K < 2 . 

The value of the variational parameter A which minimizes the latter upper bound 
may be determined from the first derivative of this expression with respect to A: 




For this particular value of A, by shuffling together all our previous inequalities, we find 
that the energy eigenvalue corresponding to the ground state of the semirelativistic 
Hamiltonian (Q) with the Coulomb potential (@), Eq, is bounded from above byEil 



En<2m\ 1-'^ 
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The reahty of this upper bound requires again k < 2. 

It is straightforward to improve this crude estimate by direct evaluation of Eq. (^. 
For our choice of trial functions, the resulting upper bound on the £Xound-state energy 
Eq may be expressed in terms of the hypergeometric function 



En < 



128 
157 




A 



Quite obviously, in this case the minimizing value of the variational parameter A must 
be determined numerically. As before, the limit A — oo tell us that now 



K < 



16 
37 



3. Numerical Approximation Methods 

For an arbitrary interaction potential, it is, in general, not possible to find a closed 
solution to the spinless Salpeter equation. Therefore, several numerical approxm^ 
methods for the (numerical) solution of this equation have been introduced.li^Ji3ll§EJ'E3 
All of these numerical schemes aim at the conversion of the spinless Salpeter equation 
into an equivalent matrix eigenvalue problem, and confine themselves to the case of a 
spherically symmetric potential V{r). Here, we briefiy sketch the maybe most efficient 
among them, namely, the "semianalytical matrix method" developed by Lucha et a/.t3 

For states of definite orbital angular momentum i described by some wave function 
ipl'x), we define the reduced radial wave function u{r) by 

^(X) = !M 37,^(^,0) ; 

r 

here yem{(^, 0) are the spherical harmonics of angular momentum i and projection m. 
For this reduced radial wave function u{r), an integral representationll3 of the spinless 
Salpeter equation may be found. Furthermore, since we consider only the case of equal 
(nonvanishing) masses m of the bound-state constituents, we may eliminate from the 
kinetic-energy part of this integral equation any dependence on the mass m by scaling 
the radial variable r like 

X := mr . 

Introducing the scaled reduced radial wave function 




as well as the dimensionless energy eigenvalue 

m 
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and the dimensionless interaction potential 

_V{r) 



V{x) :-- 



m 



m 



this particular integral representation of the spinless Salpeter equation becomesi 







E - V{x) 



U[X] 



TT 



dy Ge{x,y) 



<f £(£+1) 



+ 



+ 1 



u{y) 



(4) 



where the kernel is defined by 



with 



_ 2,2 

s = X + y , 



z = 2xy . 

Here Ki is the modified Bessel functionlll of the second kind of order £. 
In the free case, i. e., for 

V{r) = , 

as well as for potentials which are less singular than the Coulomb potential, i. e., 

V(r) oc — with 77 < 1 , 

the reduced radial wave function u{r) behaves for small r, r ^ 0, asymptotically like 

u{r) oc r^'^^ . 

Accordingly, we make the ansatz 

u{x) = x^^"^ w{x) . 

With this substitution, the integro-differential equation (11) is equivalent to 



E - V{x)] x^+^ w{x) = - dy y) y^+^ 



TT 



^2 2{i+l) d 



dy^ 



y dy 



w{y) . 



In order to solve this scaled form of the spinless Salpeter equation, we expand the 
solution 'w{x) of this integro-differential equation into a complete orthonormal system 
n = 0, 1, 2, . . .} of basis functions for L2{E?): 



N 



n=0 
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with some set of real coefficients A„. For finite A^, tliis expansion represents, of course, 
only an approximation to the exact solution w{x). The crucial pointEj of the present 
approach is our rather sophisticated choice of the basis functions which, at least in 
principle, allows for a thorough analytical treatment of the spinless Salpeter equation: 



fn{x) := V2 exp{-x) Ln{2x) , 

where Ln{x) are the Laguerre polynomials.^ 

After some straightforward algebra, the scaled spinless Salpeter equation may be 
cast into the form of a matrix eigenvalue equation for the coefficient vector A = {A„} 
in the above expansion of the solution w{x). In self-explanatory matrix notation, this 
eigenvalue equation is given bj 



A 



(5) 



where the "power matrix" and the "potential matrix" Vji^ are defined by 



p(^) 



dx X /^(x) fm{^) -^rnn 



dxx'^'V{x)U{x)Ux)=V^l 



and the "kinetic matrix" represents the action of the kinetic term in the spinless 
Salpeter equation on the vector of basis functions /„. In this way, the solution of the 
spinless Salpeter equation can be reduced to a simple matrix eigenvalue problem. The 
eigenvalues of this equation are the energies E of the bound state under consideration. 
The corresponding eigenvectors {A^} give the radial wave functions u{r) according to 



u{r) 



N 
n=0 



Kfnifnr) 



For N = 00 this treatment would be exact. For < 00 it provides an approximation 
to the exact solution of increasing accuracy with increasing N, that is, with increasing 
size of the involved matrices. As is evident from the above construction, this procedure 
works analytically, at least for all potentials of the type "power times exponential," 
i. e., for all those potentials which involve only terms of the form 



r"exp( 



with (maybe vanishing) constants 

Let us illustrate this simple prescription for the case of bound states with vanishing 
orbital angular momentum i of their constituents, that is, i = 0. All we need are the 
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explicit expressions of the kinetic matrix T^^^ and of the inverse of the nower matrix 
p(o) for i = 0. The first few entries Ti52 in the kinetic matrix T^°) readliS 



where the matrix S is given by 



TT (2n + 2m + 3)!! 



S„ 



s 



/ 1 -3 -5 -21 

-1 81 -375 -1029 

-1 -225 13125 -77175 

-3 -441 -55125 3565485 



(Note the range of the indices n and m: n,m = 0, 1, 
the power matrix P^^^ read exphcitlylll 



N.) The elements of 



oo 

1 f 

Pnl = 2]^^^ exp(-x) Ln{x) Lm{x) 


2 n + 1 for m = n 

1 I — m for m = n + 1 

— n for m = n — 1 

else 

/ 1-1 ■■■ \ 

-1 3-2 

0-2 5-3 
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V 



The inverse of this matrix, required for the matrix eigenvalue equation 
explicitly on the size N of the involved matrices and is given bytH 



depends 



forn = 0,l,...,iV , 



(p(o))-^ = V-^ 
fp(o))"^ = fp(o))"^ =(pmY^ for n = 0,1,..., AT, m = 0,l. 



n — 1 . 



For instance, for = 3 this inverse reads 



(p(o))-^ =i 

V /4x4 fi 





/ 25 


13 


7 




1 


13 


13 


7 


3 


6 


7 


7 


7 


3 




I 3 


3 


3 


3 / 
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Consider, for example, the so-called funnel potential l^(r), which depends on just 
two parameters, namely, on the Coulomb coupling constant k and on the slope a of the 
linear term: V{r) = Vp{r) = —n/r + ar. This funnel-shaped potential represents the 
prototypes of all "realistic," that is, phenomenologically acceptable, "QCD-inspired" 
static interquark potentials proposed for the description of hadrons as bound states of 
(constituent) quarks in the framework of potential models, the quarks inside a hadron 
being bound by the strong interactions arising from quantum chromodynamics. Now, 
according to the above definition, the scaled form of the funnel potential, Vy{x), reads 



which entails for the corresponding potential matrix 

In particular, specializing again to the case £ = 0, the matrix elements of the funnel 
potential, taken with respect to £ = states, are given by 

i/w = _^p(-i) + 4^p(i) . 

Since, according to the orthonormalization condition for the basis functions {fn{x)}, 
the power matrix p(~^) is identical to the unit matrix. 




p(-i) = 1 



AfxAf , 



we have 



Working out in detail the £ = 1 power matrix P^J^, we obtain! 




Explicitly, the first entries of P*^^-* read 



/ 



1 
2 
1 




2 
7 
8 
3 



1 
-8 
19 
18 




3 
18 
37 




V 



In summary, shuffling together the above results 
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• for the kinetic matrix T^^\ 

• for the inverse of the power matrix PW, and 

• for the potential matrix V^^^ of the potential under consideration 

one ends up with a well-defined matrix eigenvalue equation, which represents the 
i = special case of the general matrix form (H) of the spinless Salpeter equation. 

4. Effectively Semirelativistic Hamiltonians of Nonrelativistic Form 

Almost all of the troubles which one encounters when trying to apply the spinless 
Salpeter equation are obviously brought about by the nonlocality of the "square-root" 
operator of the relativistic kinetic energy, ^/^~+ln^, in the Hamiltonian H, Eq. (|l]). 
In contrast to the nonrelativistic limit, obtained from the expansion of the square root 
up to the lowest p^-dependent order, i/p^ + = m + p^/(2m) -|- . . ., the presence 
of this relativistic kinetic-energy operator prevents, in general, a thoroughly analytic 
discussion; one is forced to rely on numerical solutions of the problem. This (intrinsic) 
difficulty of any (semi-)relativistic formalism may be circumvented by approximating 
the Hamiltonian (||) by the corresponding "effectively semirelativistic" Hamiltonian,Bj 
which retains apparently the easier-to-handle nonrelativistic kinematics but resembles 
its relativistic counterpart to the utmost possible extent by replacing some of its basic 
parameters by effective ones which depend, in a well-defined manner, on the square of 
the relevant momentum p. The main ideaB§ of the constructionlll of these "effectively 
semirelativistic" Hamiltonians is as follows. 

Applying again the above-mentioned fundamental inequality 



which holds for any self-adjoint operator O = and arbitrary Hilbert-space vectors 
I ) (in the domain of O) normalized to unity, to the relativistic kinetic-energy operator 
y/p^~+~mF yields 



By employing this inequality, we obtain for an arbitrary expectation value (H) of the 
semirelativistic Hamiltonian H, Eq. (|I]), 



1(0)1 <v^, 





< 2 ^(p2) + m2 + {V) 
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From now on we specify the generic Hilbert-space vectors in all expectation values 
to be the eigenstates of our Hamiltonian H . In this case the expectation value of H, 
(H), as appearing, e. g., in Eq. becomes the corresponding semirelativistic energy 
eigenvalue E, i. e., 

and the inequality (|^) tells us that this energy eigenvalue is bounded from above bj 

E < (2 ^ = + V) . 



The operator within brackets on the right-hand side of this inequality may be regarded 
as some "effectively semirelativistic" Hamiltonian ifgff which possesses, quite formally, 
the structure of a nonrelativistic Hamiltonian ,!'§ 

^ (p2) + m2 

= 2m+^ + Kff, (7) 
m 

nn 

but involves, however, the effective masi 

1 



rh = (p2) + m? 



as 



well as the effective nonrelativistic potentiaflS 



2 



2 Tfi 

Kff = I -^{p^) + m? + V 

(p2) + m? 



= 2m-'^ + V . (9) 
m 

The effective mass m given by Eq. (||) and the constant, i. e., coordinate-independent, 
term in the effective potential Vcs of Eq. (|]), 

2m- ^ , 

m 

obviously depend on the expectation value of the square of the momentum p, that 
is, on (p^), and will therefore differ for different energy eigenstates. 

Motivated by our above considerations, we propose to approximate the true energy 
eigenvalues E of the semirelativistic Hamiltonian H of Eq. (|l|) by the corresponding 
"effective" enejgy eigenvalues i?efr, defined as the expectation values of some effective 
Hamiltonian iYgg taken with respect to the eigenstates |)efr of its own, 

EcS = {Hcs)cS , 



12 



where the effective Hamiltonian Hcs, as far as its structure is concerned, is given by 
Eqs. (|^ through but is imphcitly understood to involve the expectation values of 
with respect to the effective eigenstates |)efj (that is, (p^)eff in place of (p^)):Bil 

fv . ~ , - (P^)cfr , . . 

iJeflf = 4 m H z h V , 

m 

with 



^ = 2\J{V''^)eS + rii? . 

Accordingly, the effective energy eigenvalues E^s are given by a rather simple formal 
expression, viz., byBil 

Eeff = 4m+(y)eff . (10) 

We intend to elaborate our general prescription for the construction of effectively 
semirelativistic Hamiltonians Hes in more detail for the particular case of power-law 
potentials depending only on the radial coordinate r = |x|, i. e., for potentials of the 
form 

V{r) =ar'^ 

with some constant a. The reason for this restriction is twofold: 



1. On the one hand, for power-law potentials the general virial theoremci^'CJ in its 
nonrelativistic formifl appropriate for the present case, 

/p!\ 

\ ^ /eff 2 \ dr /^g ' 

enables us to replace the expectation value of the potential in ([lOD immediately 
by a well-defined function of the expectation value of the squared momentum: 

/ n\ 2 (p^)cfr 

a )efr = z — . 

n m 

This implies for the effective energy eigenvalues 

Eeff = 4m H z — . (11) 

n m 



2. On the other hand, we take advantage of the fact that for power-law potentials 
it is possible to pass, without change of the fundamental commutation relations 
between coordinate variables and their canonically conjugated momenta, from 
the dimensional phase-space variables adopted at present to new, dimensionless 
phase-space variables and to rewrite the Hamiltonian in form of a Hamiltonian 
which involves only these dimensionless phase-space variables.H The eigenvalues 
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e of this dimensionless Hamiltonian are, of course, also dimensionless.i Applying 
this procedure, we find for the effective energy eigenvalues 

m \m 



eS 



m" 



1 

2+n 



Combining both of the above expressions for E'eg, we obtain a relation which allows us 
t(|jietermine (p^)efr unambiguously in terms of the dimensionless energy eigenvalues 

= \ (^)'^"a^6^+"((p^)eff + m^) . (12) 

For a given power n this equation may be solved for (p^)efT- Insertion of the resulting 
expression into Eq. (|ll|) then yields the corresponding eigenvalue E^s of the effectively 
semirelativistic Hamiltonian H^s. 

Bound-state solutions are usually characterized by some radial quantum number 
and orbital angular-momentum quantum number £. For instance, for the harmonic 
oscillator, i. e., for n = 2, the dimensionless energy eigenvalues e are given byi 

e = 2N 

with 

3 

N = 2n, + 1+ - . 

2 

For the Coulomb potential V{r) = —n/r, that is, for n = —1, Eq. (^) reduces to a 
linear equation for the expectation value (p^)efr- Inserting the well-known expressionB 
for the dimensionless energy eigenvalues e of the (nonrelativistic) Coulomb problem. 



(2Ar)2 

with 

N = n, + e+l 
we obtain from this linear equation for (p^)efr 



2 2 

K m 



cS 



16 - ^2 ' 

and, after inserting this expression into Eq. (pH]), for the effective energy eigenvalueill 

m SN'^ - 
Epff = — 



K 



2 
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The capability of this effective formahsm to imitate the semirelativistic treatment 
becomes apparent by inspecting, for instance, the behaviour of the energy eigenvalues 
Ees for large orbital angular momenta i. For the linear potential V{r) = ar, that is, 
for n = 1, the effective formalism yields (in the ultrarelativistic limit, i. e., for m = 0) 
linear Regge trajectories:cil 

E^s = 9ai . 

This fits very well to the exact semirelativistic result ,00 
but is in clear contrast to the nonrelativistic approach, which givei 




5. Summary and Conclusion 

In spite of the regrettable fact that, at present, not even for the Coulomb potential 
the exact solution for the ground-state energy is known, the spinless Salpeter equation 
is certainly a useful tool for the semirelativistic description of bound states consisting 
of scalar bosons only as well as of the spin-averaged spectra of bound states consisting 
of fermionic constituents. 
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